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ABSTRACT 

Field theory at nonvanishing temperature beyond perturbation theory is discussed 
for the iV-component 0(iV)-symmetric scalar theory. We compute the effective po- 
tential directly in three dimensions using an exact evolution equation for an effective 
action with an infrared cutoff. A suitable truncation is solved numerically. We obtain a 
detailed quantitative picture of the scaling form of the equation of state in the vicinity 
of the critical temperature of a second order phase transition. 

1. Introduction 

The most prominent application of field theory at nonzero temperature in particle 
physics concerns the symmetry restoration in spontaneously broken gauge theories at 
high temperature. In the standard model of electroweak interactions gauge bosons 
and fermions acquire a mass due to a nonvanishing vacuum expectation value of 
the Higgs doublet. According to arguments first given by Kirzhnitz and Linde 1 this 
expectation value vanishes at sufficiently high temperature and the electroweak sym- 
metry is restored. This phase transition^ is expected to have occured in the early 
universe which may have important consequences like the possible creation of the 
excess of matter compared to antimatter (baryon asymmetry) 5 . Another challenge is 
the study of QCD at nonvanishing temperature which is also relevant for heavy ion 
collisions in the laboratory. At temperatures around a few hundred MeV one expects 
dramatic changes in the structure and symmetries of hadronic matter 6 . Of course, 
there are other fields of application of nonzero temperature field theory as for instance 
condensed matter physics. 

The prototype for investigations concerning the symmetry restoration at high 
temperature is the iV-component scalar field theory with 0(iV)-symmetry. For N = 
4 it describes the scalar sector of the electroweak standard model in the limit of 
vanishing gauge and Yukawa couplings. It is also used as an effective model for the 
chiral phase transition in QCD in the limit of two quark flavors 7 . In condensed matter 
physics N = 3 corresponds to the well known Heisenberg model used to describe the 
ferromagnetic phase transition. There are other applications like the helium superfluid 
transition (N = 2), liquid- vapor transition (N = 1) or statistical properties of long 
polymer chains (N = 0). 8 



*Talk given at the International School of Subnuclear Physics, 34th Course: Effective Theories 
and Fundamental Interactions, Erice-Sicily, 3-12 July 1996. 

**For large Higgs boson masses around m# ^ 80 GeV there might be no phase transition but 
rather an analytical crossover. 2 ' 3,4 
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The equation of state (EOS) for a magnetic system is specified by the free energy 
density per volume (here denoted by U) as a function of arbitrary magnetization <fi 
and temperature T. All thermodynamic quantities can be derived from the function 
U((f>,T). For example, the response of the system to a homogeneous magnetic field 
H follows from dU /d<p = H . This permits the computation of <fi for arbitrary H and 
T. There is a close analogy to quantum field theory at nonvanishing temperature. 
Here U corresponds to the temperature dependent effective potential as a function 
of a scalar field <fi. For instance, in the 0(4)-symmetric model for the chiral phase 
transition in two flavor QCD the meson field <fi has four components. In this picture, 
quark masses are associated with the source H ~ m q and one is interested in the 
behavior during the phase transition (or crossover) for H ^ 0. The temperature 
and source dependent meson masses and zero momentum interactions are given by 
derivatives of U. 7 

The applicability of the 0(iV)-symmetric scalar model to a wide class of very differ- 
ent physical systems in the vicinity of the critical temperature T cr at which the phase 
transition occurs is a manifestation of universality of critical phenomena. There exists 
a universal scaling form of the EOS in the vicinity of a second order phase transition 
which is not accessible to standard perturbation theory. The quantitative description 
of the scaling form of the EOS will be the main topic of this talk. To compute the 
effective potential U which encodes the EOS we employ a nonperturbative method. 
It is based on an exact evolution equation for an effective action r& with an infrared 
cutoff ~ k. From the standard effective action T, i.e. the generating functional of 
IP I Green functions, is recovered when the additional infrared cutoff ~ k is removed. 
Though I will refrain from going into technical details, an introduction to the method 
and some results for the EOS 9 will be presented. 

2. Effective three dimensional behavior of high temperature field theory 

Quantum field theory at nonzero temperature T can be formulated in terms of an 
Euclidean functional integral where the 'time' dimension is compactified on a torus 
with radius T -1 . 10 If the characteristic length scale of the considered physical problem 
is much larger than the inverse temperature the compactified dimension cannot be 
resolved ('dimensional reduction' 11 ). One therefore observes an effective three dimen- 
sional behavior of the high temperature quantum field theory. If the phase transition 
is second order the correlation length becomes infinite and the effective three dimen- 
sional quantum field theory is dominated by classical statistics. In particular, the 
critical exponents which describe the singular behavior of various quantities near the 
phase transition are those of the corresponding classical statistical system. 12 Analo- 
gous considerations are valid for sufficiently weak first order phase transitions. 

The calculation of the effective potential U (0, T) in the vicinity of the critical 
temperature of a second order phase transition is an old problem. One can prove 
through a general renormalization group analysis 13 the Widom scaling form 14 of the 
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Only the limiting cases 0^0 and — > oo are quantitatively well described by critical 
exponents and amplitudes. The critical exponents (3 and 5 have been computed with 
high accuracy 8,15 but the scaling function / is more difficult to access. A particular 
difficulty for a direct computation in three dimensions arises from the existence of 
massless Goldstone modes in the phase with spontaneous symmetry breaking for mod- 
els with continuous symmetry (N > 1). They introduce severe infrared divergencies 16 
within perturbative expansions. A related example for the breakdown of perturba- 
tion theory is the symmetric phase of the electroweak theory where the gauge bosons 
are massless. In addition, within a perturbative description one is restricted to small 
expansion parameters. For instance, in the three dimensional 0(iV)-symmetric model 
no small coupling characterizes the interactions at the phase transition. One therefore 
has to employ nonperturbative methods which resolve the infrared problem. 

A standard method to cope with infrared singularities is the 'e-expansion' which 
is a double series expansion in powers of the coupling constant g and e = 4 — d. It 
has been used to compute the scaling function / in second order in e (third order 
for iV = l) 8 . To study the three dimensional theory one sets the expansion param- 
eter e — 1. The validity of this expansion is not clear a priori since the expansion 
parameter is not small. A nonperturbative computational method is provided by lat- 
tice Monte-Carlo simulations 19 . Nevertheless second order or weak first order phase 
transitions are very demanding due to a finite lattice size and an infinite or large 
correlation length. In particular, it is notoriously difficult to distinguish by these 
methods between a weak first order and a second order nature of the transition. The 
nonperturbative method which allows here an unambiguous answer and which we will 
consider in the following relies on the use of an 'exact renormalization group equa- 
tion' (ERGE). ERGEs have been formulated in many different but formally equivalent 
ways 13,17 . Though these equations are exact there is practically no chance to solve 
them without approximations. The main challenge is to find a suitable approximation 
scheme that describes all the relevant physics of the considered problem. It is this 
point where the particular formulation of an ERGE becomes important. Here we will 
use the effective average action method 18 which allows nonperturbative calculations 
in a feasible way. Before turning to the method in section 4 we consider some of its 
results for the scaling EOS for 0(iV)-symmetric models. 9 They have been obtained 
from an approximate numerical solution of the effective potential U directly in three 
dimensions (cf. sect. 4). 

3. Scaling equation of state 

Eq. (jl]) establishes the scaling properties of the EOS. They have been explicitly 
verified by our numerical solution 9 . In particular, the function / depends on only one 
scaling variable x. In our conventions x = —Ska/^ 1 ^ where Ska is proportional to 
the deviation from the critical temperature Sk\ = A(T)(T cr — T) with A(T cr ) > 
and = v '0 a a (a — 1, ... , N). If one plots the logarithm of / as a function of the 
logarithm of its argument x one can easily consider some limits and compare with 
known values of critical exponents and amplitudes. Figs. 1 and 2 show our results 
for log(/) and \og(df /dx) as a function of log|x| for N — 1 and N = 3. Here x > 
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Figure 1: Logarithmic plot of the scaling function / in the symmetric phase. 



corresponds to temperatures above the critical temperature (symmetric phase) and 
x < accordingly to T < T cr (phase with spontaneous symmetry breaking). 

As an example we consider the limit x — ► oo. One observes that log(/) becomes 
a linear function of log(x) with constant slope 7. In this limit the universal function 
takes the form 

lim f(x) = (C+)- 1 x 7 . (2) 



The amplitude C + and the critical exponent 7 characterize the behavior of the 'un- 
renormalized' squared mass or inverse susceptibility 
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We have verified the scaling relation 7//? = 8 — 1 that connects 7 with the exponents 
P and 5 appearing in the Widom scaling form (]!]). For the singular behavior of the 
susceptibility for 5ka — ► in the symmetric phase one finds x — C + \5ka\~~ 1 with 
7 = 1.258(1.465) for iV = 1(3). In another limit, x — > 0, the scaling function becomes 
a constant, /(0) = D, and according to eq. (|T|) one finds H = D<p & (5 = 4.75(4.78) 
for N = 1(3)). 

The spontaneously broken phase is characterized by a nonzero value 0o of the 
minimum of the effective potential U with H = (dU/d<p)(<j)o) = 0. The appearance of 
spontaneous symmetry breaking below T CT implies that f(x) has a zero x = —B~ 1 ^ 
and one observes a singularity of the logarithmic plot in fig. 2. In particular, according 
to eq. (|) the minimum behaves as O = B(5k k Y ((3 = 0.336(0.388) for N = 1(3)). 
The critical exponents like 7, (3 and 5 are universal whereas the scaling function 
/ is universal up to arbitrary normalizations of x and f(x). There are only two 
independent scales in the vicinity of the transition point which can be related to the 
deviation from the critical temperature and to the external source. All models in the 
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Figure 2: Logarithmic plot of the scaling function / in the spontaneously broken phase. 

same universality class can be related by a multiplicative rescaling of <fi and 5k\ or 
(T cr — T). Accordingly there are only two independent amplitudes and exponents 
respectively. 

Apart from the given examples we have extracted other quantities which char- 
acterize the asymptotic behavior for the sake of comparison. The critical exponents 
and amplitude ratios we have calculated typically deviate by a few per cent from the 
values obtained by other methods 8 ' 15 . We expect the error to be related to the size of 
the anomalous dimension rj ~ 4%. For recent lattice results for the EOS and a com- 
parison see ref. 19. A more detailed discussion of the EOS, including semi-analytical 
expressions for the scaling function f(x) for iV = 1 and N = 3 and an alternative 
parametrization of the EOS in terms of renormalized quantities, can be found in ref. 9. 

4. Effective average action and exact evolution equation 

In the following we turn to the method we have used to derive the above results. 
It relies on the effective average action 18 T k . For a theory described by a classical 
action S, the effective average action results from the integration of fluctuations with 
characteristic momenta larger than a given infrared cutoff ~ k. To be explicit we 
consider a /c-dependent generating functional with N real scalar fields Xa, 

W k [J] = In/ D X exp (-S\x] - A k S[ X ] + / d d xJ a {x) X a {x)^ (4) 

where the k-dependence arises from the additional term 

AkS[x] = \ I ^- d R ^)Xa(-q)x a (q)- (5) 

Without this term W k becomes the usual generating functional for the connected 
Green functions. Here the infrared cutoff function R k is required to vanish for k — > 
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and to diverge for k — > oo and fixed q 2 . This can be achieved, for example, by 
the choice Rk(q) = Z k q 2 e~ q l k (1 — e~ 9 l k ) _1 where Z k denotes an appropriate wave 
function renormalization constant which will be defined below. For fluctuations with 
small momenta q 2 <C k 2 the cutoff R k — Z k k 2 acts like an additional mass term and 
prevents their propagation. For q 2 ^> k 2 the infrared cutoff vanishes such that the 
functional integration of the high momentum modes is not disturbed. The expectation 
value of x m the presence of A k S[x] and J reads 4> a =< \ a >— $W k [J]/5J a . We 
define the effective average action as 

r fe [0] = -W k [J] + J d d xJ a {x)4> a {x) - A k S[<j>]. (6) 

As a consequence, as the scale k is lowered T k interpolates from the classical action 
S = lim^ooFfc to the standard effective action T = lim^or*;, i.e. the generating 
functional of IP I Green functions. 18 Lowering k results in a successive inclusion of 
fluctuations with momenta q 2 ^k 2 and therefore permits to explore the theory on 
larger and larger length scales. One can view T k as the effective action for averages 
of fields over a volume of size ~ k~ d and the approach is similar in spirit to the 
block spin action 20 ' 13 on the lattice. The interpolation property of T k can be used 
to 'start' at some high momentum scale A where T\ can be taken as the classical or 
short distance action and to solve the theory by following T k to k — > 0. The scale 
dependence of T k can be described by an exact evolution equation 18 

!r, M = ia>{(i? M + *)-'^} (T) 

where t = ln(/c/A). The evolution is described in terms of the exact inverse propagator 
(r^)%(q,q') = 5 2 r k /5<f)a{—q)84> b (q'). The trace involves a momentum integration as 
well as a summation over internal indices. The additional cutoff function R k with a 
form like the one given above renders the momentum integration both infrared (IR) 
and ultraviolet (UV) finite. In particular, the direct implementation of the additional 
mass term R k ~ Z k k 2 for q 2 <C k 2 into the inverse average propagator makes the 
formulation suitable for dealing with theories which are plagued by infrared problems 
in perturbation theory. 

Though the evolution equation for the effective average action is exact it remains a 
complicated functional differential equation. In practice one has to find a truncation 
for T k and obtains approximate solutions. Our ansatz represents the lowest order in 
a systematic derivative expansion of T k , 

r k = j d d x {U k (p) + z k d^ a dy a } , a — 1, ... ,N. (8) 

Here (f> a denotes the A^-component real scalar field and p = \4> a 4> a - We keep for the po- 
tential term the most general 0(A r )-symmetric form U k (p) since U(p) = lim fc ^ U k (p) 
encodes the equation of state. The wavefunction renormalization is approximated 
by one /c-dependent parameter Z k . The first correction to this ansatz would include 
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field dependent wave function renormalizations Zk(p) plus functions not specified in 
eq. (HD which account for a different index structure of invariants with two deriva- 
tives for N > l. 12 ' 18 The next level involves invariants with four derivatives and so 
on. Concerning the equation of state the truncation of the higher derivative terms 
typically generates an uncertainty of the order of the anomalous dimension 77. 9 For 
N = 1 the weak p-dependence of Z k has been established explicitly at the critical 
temperature 21 . We define Z k at the minimum p of U k and at vanishing momenta q 2 , 
i.e. Z k = Z k (p = p ; q 2 = 0). The fc-dependence of this function is determined by the 
anomalous dimension 7](k) = — d(lnZ k )/dt. 

If the ansatz is inserted into the evolution equation for the effective average 
action (|7j) one obtains flow equations for the effective average potential U k (p) and 
for the wave function renormalization constant Z k (or equivalently the anomalous 
dimension rj). These have to be integrated starting from some short distance scale A 
and one has to specify U\ and Z\ as initial conditions. The short distance potential 
is taken to be a quartic potential 

U K {p) = -p 2 k p+ l -\ KP 2 (9) 

and Z\ = 1. For a study of the behavior in the vicinity of the phase transition, it is 
convenient to work with dimensionless renormalized fields p = Z k k 2 ~ d p and we also 
switch to a dimensionless potential u k = k~ d U k . The evolution equation for u k (p) 
reads 12,18,21 (primes denote derivatives with respect to p at fixed t) 

= - duk + (d - 2 + r})pu' k + 2v d (N - l)l d (u' k - rj) + 2v d l d (u' k + 2p<; V ) (10) 

where v^ 1 = 2 d+1 ir d ^ 2 T(d/2) with 1)3 = 1/87T 2 . The anomalous dimension rj is given 
in our truncation by 12,18 

V (k) = ^ K «(k))X, 2 (2«<(k) ;) )) (11) 

with k the location of the minimum of the potential, u' k (n) = 0. The functions l d 
and m d 2 result from the momentum integration in the flow equation (0) for T^. They 
equal constants of order one for vanishing arguments and decay fast for arguments 
much larger than one. 12,22 As a consequence these functions account for the decou- 
pling of modes with masses much larger than the scale k. For the minimum of the 
dimensionless potential at k 7^ one easily recognizes in eq. (|l~0|) the contribution due 
to the (N — 1) massless Goldstone bosons with dimensionless mass term u' k (n) = 
and the radial mode with dimensionless mass term 2ku'1(k). 

At a second order phase transition the correlation length diverges and accordingly 
there is no mass scale present in the theory. In particular, one expects a scaling 
behavior of the effective average potential. Exactly at a second order phase transition 
u k should be given by a ^-independent (scaling) solution du k /dt = 0. The EOS 
involves the potential away from the critical temperature. Its computation therefore 
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Figure 3: The evolution of u'Ap) as k is lowered from A to zero. 

requires the solution of the full partial differential equation (|T0|) for the dependence 
of Uk on the two variables p and k. Suitable numerical methods for its solution are 
discussed in ref. 22. 

In the remaining part of this section we consider the theory in three dimensions 
(d = 3) with the classical potential ([D as initial condition. For given Aa/A there 
is a critical value for the minimum n\ = p\/(X\A) = n cr of the classical potential 
for which the scaling solution is approached in the limit k — > 0. In particular, the 
minimum n and the anomalous dimension rj take on constant (fixed point) values 
n(k) = k ± and rj(k) = rj*. The order of the phase transition becomes apparent from 
the observation that the order parameter po ~ hrcifc_*o k 1+v *K+ = vanishes at the 
phase transition. 

To demonstrate the behavior of the scale dependent effective potential near the 
phase transition we consider the shape of u' k (p) for different values of the scale. In 
fig. 3 the numerical solution of the function u' k (p) is plotted for d = 3 and N = 1 
and various values of t = ln(fc/A). 22 The evolution starts at k = A (t = 0) where 
no integration of fluctuations has been performed. We arbitrarily choose Aa/A = 0.1 
and fine tune that the scaling solution is approached at later stages of 

the evolution. If k\ is interpreted as a function of temperature the small deviation 
proportional to the deviation from the critical temperature. For 
the example in fig. 3 a value slightly above the critical temperature is used (<5ka < 0). 
As k is lowered u' k (p) deviates from the initial linear shape of u' A (p). By staying near 
the scaling solution for several orders of magnitude in k, the system looses mem- 
ory of the initial conditions at the short distance scale A. As a result, after u' k (p) 
has evolved away from the scaling solution, its shape is independent of the choice of 
Aa/A for the classical theory. This property gives rise to the universal behavior near 
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second order phase transitions discussed in section 3. Eventually, for the final part 
of the evolution as k — > 0, the theory settles down either in the symmetric phase 
(Ska < 0) as is demonstrated in fig. 3 with k — and positive constant mass term 
m 2 = linifc^o k 2 u' k (0), or k grows in such a way that po = linifc^o Z k x kK approaches 
a constant value indicating spontaneous symmetry breaking (Ska > 0). The scaling 
EOS obtains from U'(p) = \im k ^ Z k k 2 u' k (p) with H/(p s = U'(j) 5 ~ l = f as a function 
of x = -Ska/<P 1/p (0 = yflp). 

5. Conclusions and outlook 

We have arrived at a detailed quantitative picture of the equation of state in the 
vicinity of the critical temperature of a second order phase transition. The method we 
have used is free of IR or UV divergencies and provides a practical tool for calculations 
not accessible to perturbation theory. It relies on the effective average action F k which 
interpolates between the classical action S (k — A) and the effective action T (k = 0). 
Its scale dependence is described by an exact nonperturbative evolution equation. 

Though not discussed in the present talk, the method can be applied as well 
to a description of first order phase transitions along the lines presented above. In 
particular, the scaling equation of state for weak (fluctuation induced) first order 
phase transitions in scalar matrix models has been obtained. 23 Work in progress 
focuses on the equation of state for an effective quark-meson model 24 at nonzero 
temperature 25 , relevant for the study of the chiral phase transition in QCD. 

The average action approach allows a well-defined description of a rich spectrum 
of critical phenomena like the observation of critical exponents, tricritical points and 
crossover behavior. Further applications at nonzero temperature in various levels of 
truncations include studies of the crossover behavior in two-scalar theories 26 or the 
description of the Kosterlitz-Thouless phase transition in two space dimensions 27 . 
Extensions to gauge theories 3 ' 28 ' 29 have been e.g. applied to the abelian Higgs model 
and to the electroweak phase transition at nonzero temperature 3,28 . As should be 
pointed out the method is not restricted to a discussion of the universal behavior of 
the theory which can be observed near or at second order or sufficiently weak first 
order phase transitions. It also permits the description of nonuniversal physics such 
as the behavior of the theory away from the critical temperature for the four dimen- 
sional theory at nonvanishing temperature 12 . 
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